IAES International Journal of Robotics and Automation (IJRA) 
Vol. 11, No. 2, june 2022, pp. 132~140 
ISSN: 2722-2586, DOI: 10.1159 1/ijra.v11i12.pp132-140 g 132 


Kinematic analysis and design of a new two-limb parallel 
Schönflies-motion generator considering isotropic 
configuration 


Guanglei Wu!, Huiping Shen? 
1School of Mechanical Engineering, Dalian University of Technology, Dalian, China 
?School of Mechanical Engineering, Changzhou University, Jiangsu, China 


Article Info ABSTRACT 
Article history: This paper presents a new two-limb parallel Schénflies-motion generator, which 
Received Feb 10. 2021 adopts a pair of alternative spatial modules equivalent to the parallelogram struc- 


ture. This modular architecture can ensure the enhanced stiffness of the manipu- 
lator normal to the motion of the planar parallelogram structure due to the trape- 
zoidal architecture. The preliminary kinematic problems, namely, the mobility, 
forward/inverse geometry and singularity, are studied as well as the kinematic 
Keywords: isotropy. A neutral isotropic configuration of the robot is identified for the struc- 
tural design of the link lengths. 


Revised May 25, 2021 
Accepted Jun 23, 2021 


Kinematic isotropy 


Parallel robot 

Schénflies-motion generator This is an open access article under the CC BY-SA license. 

= 
BY SA 

Corresponding Author: 

Guanglei Wu 


School of Mechanical Engineering, Dalian University of Technology 
116024 Dalian, Liaoning, China 
Email: gwu @dlut.edu.cn 


1. INTRODUCTION 

Parallel robots with Schönflies-motion (three independent translations and one rotation around a 
fixed axis), also known as, selective compliance assembly robot arm (SCARA) motion [1], named as parallel 
Schönflies-motion generators, are more and more attractive for high-speed pick-and-place (PnP) applications 
serving light industries, thanks to their lightweight architecture and high stiffness. Amongst this type of robots 
[2]-[8], most of them so far have the symmetric structure with four identical limbs, which admit a cylinderical 
workspace. On the other hand, the PnP operations of the robots deployed in the production line are usually 
confined in a rectangular region. In view of the properties of PnP motion, the PnP robots with rectangular 
workspace are desired, allowing more robots to be installed side-by-side, to utilize the shop-floor space more 
efficiently. For this purpose, besides the recently developed asymmetrical four-limb robot [9], the architecture 
of the two-limb Schönflies-motion generators (SMGs) can effectively reduce the mounting space. 

The design and development of the two-limb SMGs have been reported in the literature. For instance, 
Angeles et al. [10] developed the McGill SMG with two parallelogram structures (II joint) in one limb, which 
has a rectangular mounting space. On the other hand, the second floating motor in each limb leads to the 
increased inertia. Later, Lee and Lee [11] proposed several asymmetric SMGs with fewer links, which behaves 
isotropic kinematic and static performances. Harada and Angeles [12] selected one robot architecture from the 
previous work for prototyping, named as “Pepper Mill Carrier”. Although this robot architecture simplified 
the robot linkage design, the complexity of the actuation system increases due to the left/right-hand ball-screw 
driven structure. Some other designs inspired from this concept have also been reported [13], [14]. It should 
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be noted that this asymmetric architecture results in the property of unbalanced dynamics [15]. Another design 
of the five-bar based two-limb SMG [16] has fixed motor to reduce the inertia, but with decreased stiffness and 
limited motion along vertical axis. Linear actuated two-limb SMG developed by Wu et al. [17] can have almost 
the same workspace along its linear rail of the actuator. The previous two-limb SMGs saved shop-floor space 
due to their architecture, whereas, the stiffness along some motion axis decreased, unlike the fully symmetrical 
four-limb SMGs. Upon the investigation of the stiffness behavior [18], the stiffness along the normal to the 
motion of the planar parallelogram structure can be improved with an alternative spatial module, which can 
have different architectures [19], [20], towards high dynamic response. Table 1 lists a brief comparison of the 
reported two-limb SMGs. This paper will present the kinematic analysis of a new two-limb parallel Schönflies- 
motion generator based on the lightweight design of the pair of spatial modules. 

In this paper, a new two-limb parallel Schönflies-motion generator is introduced, which adopts a pair 
of alternative spatial modules equivalent to the parallelogram structure. This modular architecture can ensure 
the enhanced stiffness of the manipulator normal to the motion of the planar parallelogram structure due to 
the trapezoidal architecture. The kinematic problems, namely, the mobility, forward/inverse geometry and 
singularity, are studied as well as the kinematic isotropy. A neutral isotropic configuration of the robot is 
identified for the structural design of the link lengths. 


Table 1. Comparison of the reported two-limb SMGs 


SMGs Advantages Disadvantages 
McGill SMG [10] Fewer links Increased dynamic inertia due to movable mo- 
tors 
C-drived SMGs [12]-[14] Fewer links, isotropic kinematic and Inherently unbalanced dynamics, complex cylin- 
static performances drical joint design 
Five-bar based SMG [16] Fewer links, low dynamic inertia Decreased stiffness, limited motion along verti- 
cal axis 
Linear SMG [17] Large workspace along its linear rail of Decreased stiffness along the linear rail 
the actuator 
Proposed SMG Enhanced stiffness, fewer linkages Trapezoidal lower linkages, more kinematic 
joints 


2. RESEARCH METHOD 
2.1. Manipulator architecture 

Figure 1 depicts the conceptual design of the two-limb parallel Schönflies-motion manipulator, for 
which each limb of the robot is actuated by two perpendicular actuators on the base. A pair of spatial modules, 
instead of the conventional II joints, are adopted to the connect the mobile platform and the proximal links. 

The parameterized kinematic structure of the robot is illustrated in Figure 2, where the global co- 
ordinate frame F, is built with the origin located at the geometric center of the base platform, the x-axis 
being collinear to segment A1 A2. The moving coordinate frame F, is attached to the mobile platform and 
the origin is at the geometric center, where X-axis is parallel to segment CC 2. Here and after, vectors i, 
j and k represent the unit vectors of x-, y- and z-axis, respectively. The axis of rotation of the ijth actu- 
ated joint of angular displacement 0;; E€ (—1/2, 7/2), i = 1, 2, j = 1, 2, is parallel to u;; = R(z, a;;)j, 
aij = (—1)'*'1/2 + (—1)47/4. Moreover, as the offset distance e on the elbow bracket does not affect the 
output variable, but is added as an offset to z-coordinate, henceforth, the offset distance e is assumed to be zero 
for convenience. 


2.2. Degree-of-freedom and displacement group 

The degree-of-freedom (DOF) of the robot under study is derived by the group theory, consequently, 
according to Figure 3, the bond £; of the ith limb is the product of the following bonds i) The rotation subgroup 
R(A;,;) of axis of rotation A;; passing through A;, and parallel to u;;; ii) The translation subgroup 7 (n;;) 
corresponding to the II-joint lying in a plane ¢;; normal to n;;; iii) The motion subgroup of the spatial module 
being the product of rotation subgroup R(;) of axis of rotation B; passing through B; and translation subgroup 
T (y) normal to y-axis; and iv) The rotation subgroup R(C;) axis of rotation C; passing through C; and parallel 
to k. Thus, the kinematic bonds of the ith limb is expressed as (1) and (2). 
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& = T (n;) - R(B:) -T (y) -R(Ci) = R(B:) : X (k) (1) 
Where n; is perpendicular to v;. Subsequently, the intersection of the bonds £ and £% yields. 


Lı N L2 = X(k) (2) 


Henceforth, the intersection of all subgroups being a Schénflies subgroup Æ (k), the robot generates the 
Schönflies motion. 


Elbow bracket __ 
of Limb 7 | - 


Universal 
joints 


Spatial module 
of Limb 7 


Wrist bracket of Limb / Mobile platform 


Figure 1. The conceptual design of the new two-limb parallel Sch6nflies-motion robot 


Figure 3. The loop structure of the robot with rotational inputs 
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2.3. Kinematic geometry of the manipulator 


Under the coordinate systems of Figure 2, the position vectors of point A; in frame F, are denoted by 
a; = (—1)‘ai, then the Cartesian coordinates of point B; is derived as (3). 


b; = bv; +a;; vi = Rik = Ra Riek (3) 


Where R;; = R(z, ai;)R(y, 9:;) and 6;; is the input angle from the reference plane. 
Let the mobile platform pose be denoted by x = [p7 o] D p= |z y z| p the Cartesian coordi- 
nates of point C; is expressed as (4). 
ci = Qc; +p (4) 


Where Q = R(z, ¢) is the rotation matrix of the MP and c’, = (—1)Żri is the position vector of point C; in the 
frame Fp. 


2.3.1. Forward geometry 


Due to translational motion of the distal spatial modules, the y-coordinate of point B; is equal to that 
of C;, thus, the y-coordinate of point P and the orientation of the mobile platform, respectively, are derived as 


(5). 


1 (b2 — bij 


2r 2 


2 
1 
y= 22 bij: ọ = sin" 
On the other hand, the x- and z-coordinates of point P of the mobile platform can be solved from the 
four-bar linkage as shown in Figure 4, namely, expressed as (6). 
(sing, — sing.) + gı =0; I(cosq, + cos q2) + g2 = 0 (6) 


Where gı = (bz — b1)fk and gz = 2r cos ¢ + (bı — bz)? i. Consequently, the angles q1, q2 € [0, 7/2] are 
obtained as (7). 


1 291 + VAF? + 492 — (9? + 93)? 


es We 
; q2 =sin™ (sing: + g1) (7) 
gi +93 — 292 


qı = 2tan— 


To this end, one obtains the x- and z-coordinates of point P in the frame F, is expressed as (8). 
n= cfi; z=10 ck (8) 
24 2 4 í 


, stands for the position vector of point C;. 


J7 


Where c; = b; + lw;, w; = |(—1)}+t cosq; 0 sing 


Figure 4. The planar four-bar linkage B,C,C Bo in the forward geometry 
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2.3.2. Inverse geometry 


The inverse geometric problem is solved from the following kinematic constraints expressed as (9) 


and (10). 
(ci — by) -j = 0, i=1,2 


(ci = b;)" (c; = b;) = ? 
Expanding the above equations leads to (11a) and (11b). 
Di, sin Oi + Di, cos 6; + di = 0 


D5, sin 6,1 + Dg cos 6:1 + dy = 0 


Which can be cast in a matrix form as (12). 
Be a ea _p. | on H 
Dı D3ə| |cos Ai ? [cos 641 d; 
With the coefficients expressed as (13a), (13b), and (14a) to (14e). 
Diy = hachi; Dis = hizsbi2; Dd, = hizsðiz + hi4chi2; Dho = hissbi2 + higcOj2 


di = dashi + diz; d3 = digsi2 + dia 


hii = —bsa;; hig = —bcai2SQi1 


hi3 = 2bzcaiz; hi4 = —2b( (x + a)cair + YSQji1 re(ai ¢)) 
his = —2beaj2((x + a)ca;ı + ysa; — refa; — b)); hie = —2bz 
dj, = —bcaj1saj2; dig = yY — rsp 


di3 = 2bsa;2((x + a)sair YcAil rs(&ir @)); dia = lle; = a; ||? + b? an P 


(9) 
(10) 


(11a) 
(11b) 


(12) 


(13a) 
(13b) 
(14a) 
(14b) 
(14c) 
(14d) 
(14e) 


Where s and c stand for the sine and cosine functions, respectively. The determinant of matrix D; is calculated 


as (15). 
|[D;| = 2b?°z sin a; (1 — sin? a; sin? 0;2) 4 0 


(15) 


Solving (12) using Cramer’s rule, sin 0;ı and cos 0;1 can be obtained as the functions of 6;2. Expanding 


and simplifying equation sin” 6;, + cos? 6;; = 1 yields expressed as (16) with (17a) to (17g). 
Sy sinf Big + So sin? i2 + S3 cos? Oiz sin? Oiz + S4 sin? 652 
+55 cos? Oiz sin O52 + S6 cos* 6:2 + S7 cos” Oin = 0 
Sı = & (hig + hig) — 2didizshizdis + hi (diz — his) 
S2 = 2dirdia (hig + his) + 2dizdishiz — Wi2dis(dirdia + dizdi3) 
S3 — d (h3 + his) — 2d; dighi1 dia + hey (d = his) + Qhithia(hiahis + hishie) = heh, 
Sa = dig(hig + his) — Qdiadiahiodis + hindig 
Ss = 2dirdi2 (hi, + hig) + 2digdiah?, — 2hirdia(dirdia + dizdis) 
So = -h his 
S7 = din(hig + hig) — 2diadiahirdia + djhi, 


Equation (16) can be written in the following equivalent form (18) with (19). 


sinf Oiz + Li sin? Oiz + Mi; sin? iz + N; sin Oiz + P; =0 


So — Ss eth S4 = 2S6 — S7. i Ss . p= Se + S7 
S1 — S3 + S6 i i S1 — S3 + S6’ g S1 — S3 + S6 


besoa a M= 
S1 — S3 + S6 
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As a result, angle 0;2 can be obtained by solving the Quartic (18), namely, expressed as (20) with (21), 
where expressed as (22a) to (22c). 


Li die 2M; 1 /L? 4M; 
0; — aint ? y 4 4 I; + Jj =e y 2 Li yE Ki 2 
2 = Sin ( 1 5Va4 3 + J, zV? 3 J, (20) 
o VU, pagam. = 1 -I +41;M; —8N; (21) 
t 39W V 54? 4/0274 2M, /3 41, + Si 
U; = M? — 3L,N; + 12P, (22a) 
V; = 2M? — 9L;M;N; + 27N? + 27L?P, — 72M;P, (22b) 


W; = V; + \/—4U} + V? (22c) 


Substituting ;2 into (9), angle 0;1 is then solved as (23) with (24a) to (24c). 


—X,+ X24 Y2 — Z? 
i ten või tř 4i (23) 


Zi- Yi 
Xi = —b sin aj, cos 652 (24a) 
Y; = —b cos ajg sin ayy sin ĝ;2 (24b) 
Zi = y — rsin ġ — bsin qiz cos aj sin 642 (24c) 


3. RESULTS AND DISCUSSION 
3.1. Jacobian analysis 

The Jacobian matrices of the robot can be obtained by differentiating (9) and (10) with respect to 
time. As the y-coordinate is independent of link length B;C;,, here, in order to make these two equations 
homogeneously dimensionless, (9) and (10) can be rewritten as (25) and (26). 


(ci —b;)-j/b=0, i=1,2 (25) 
(ci — bi)” (ci — b;)/? -1 = 0 (26) 
Upon the differentiations of (25) and (26), one obtains expressed as (27) and (28). 
(dr cos qi sing + wT p)/l = (ĝbw? (ua x vi) + ĝizbw? (uiz x v:))/l (27) 
((-1)'ér cos b +j” p)/b = bij” (uit x vi) + Bini" (u2 x vi) (28) 


Where p = [i y 3] ie Equations and can be cast in a matrix form, namely, expressed as (29) with 
(30a) to (30c). 
Ax = B9 (29) 


cos q1 0 sing, rcosq sing 
A= 0 l/b 0 =rl/bcosġ (30a) 
—cosqg 0 sing. rcosq sing 


0 l/b 0 rl/b cos ġ 


bw? (u11 x vı) bwT (u12 x vı) 0 0 
= Bı 02 = Ij? (u11 x vı) Ij? (ure x vı) 0 0 
B E P A T 0 0 bw (u21 X v2) bwi (u22 X v2) (30b) 
0 0 Ij? (u21 x V2) Ij? (u22 x v2) 
. à 5 á -T : : r š . -T 
x=[¢ 9 2 d) ; ð=[ġ ð> ðs 64] (30c) 


Where A and B are the forward and backward Jacobians, respectively. As along as A is nonsingular, the 
kinematic Jacobian matrix is obtained as (31). 


J = A`!B (31) 
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3.1.1. Singularity analysis 

The three types of singularities can be determined from the Jacobian matrices in (29) namely, i) 
det(B) = 0, the robot loses one or more DOF and reaches a serial singularity; ii) det(A) = 0, the robot gains 
one or more uncontrolled DOF and encounters a parallel singularity; and iii) det(A) = det(B) = 0, the robot 
is in a configuration of mixed singularity. Here, we only focus on the first two types of singularities, as this 
type of singularities can be obtained from the two previous ones. 

From det(B) = det(B;) - det(B2), the robot reaches type 1 singularity in the configuration uj; || v; 
or/and w; || v;, which usually occurs at the workspace boundaries. The type 2 singularities arise when w; || 
Wo, Le. q1 = q2 = 0 or qı = q = 7/2 or 6 = +r /2, which leads to the uncontrollable motion of the mobile 
platform. 


3.1.2. Isotropic configuration 
The manipulator isotropy is often related to the dexterity, and the Jacobian matrix is said to be isotropic 
under the conditions expressed as (32). 


ATA=7°L; B'B = 0°, (32) 


Where 7 and o denote the four identical singular values that each Jacobian matrix must have at the optimum 
postures, and I in the identity matrix. As the entries of the forward Jacobian matrix of the robot are not 
homogeneous as they contain mixed rotation and translation terms, for which the condition number is limited 
to indicate the kinematic performance. Here, the characteristic length is introduced to normalize the forward 
Jacobian matrix [21], where the characteristic length Le is obtained as (33). 


Ax=[A, A] A ={A, A,/Le| | 


Where A, is the first three columns and A, is the last one in A, respectively, and A is the dimensional 
homogeneous Jacobian. Consequently, the isotropy condition can be expressed as (34). 


ATA, Rr N E °] 


D | a 
P| = AX (33) 


ATA, /Le ATA,/L2|~ | 0? 7 G4) 


T 


ATA=| 


Hence, the characteristic length is obtained as (35). 


3ATA 
L.= ew 
\ tr(ATA,) Ga) 


Substituting (35) into (34) yields expressed as (36) with the characteristic length and the identical singular 
value, respectively, expressed as (37). 


q =Q =7T/4 (36) 
LH Ts THT (37) 


The backward Jacobian matrix B is isotropic when BTB = o71y, consequently, the common singular 
value is found as (38) with the isotropic configurations as (39). 


{ o=1/V5; a =0, —O12 = 622 = tan~12 (38) 


et ee. g=1 (39) 


g=0,y¥=0, 2=b/V541/V2,¢=0; o =1/V5 

In accordance with (39), it is seen that the isotropic configurations are related to both manipulator 

pose and the geometric parameters. For the robot under study, its isotropic configurations exist when the 
design parameters meet the following conditions as (40). 


l= V2(a—r),b=V2l or 1=V2(a—r—b/V5), b = V2 (40) 
In practice, the first case, i.e., 1 = V2(a —r),b= J/2l, can be selected as the design candidate to define the 


home configuration of the robot because of 6;; = 0, which is in an isotropic configuration. 
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4. CONCLUSION 

In this paper, a new two-limb parallel Sch6nflies-motion generator is introduced, which adopts a pair 
of alternative spatial module equivalent to the parallelogram structure. Compared to the fully symmetrical 
four-limb counterparts, this robot architecture can help to save the shop-floor space, allowing more robots to be 
deployed in the production line side-by-side, to improve the productivity. Moreover, this modular architecture 
can ensure the enhanced stiffness of the manipulator normal to the motion of the planar parallelogram structure. 

The mobility of the robot is verified with the Lie group, which shows that the robot can generate 
Schonflies-motion. Similar to other types of two-limb SMGs, the forward geometry problem of the new con- 
figuration is simple, and the inverse geometric problem can also be solved symbolically. The manipulator sin- 
gularity is studied by means of the Jacobian method, where singular configurations are identified. Moreover, 
the kinematic isotropy is investigated in pursuit of the identical singular values of the homogenous Jacobian 
matrices. An isotropic configuration of the robot, which can be considered as the neutral configuration, is de- 
tected, showing the geometric relationship of the link lengths for the robot structural design. The mechanical 
design of the robot prototype is under way for experimental validation and testing. Future work will focus on 
the robot dynamics and control. 
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